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Abstract 

We have studied the magnetophoretic particle separation and its nonUnear behavior of ferrofluids in mi- 
crochannel which is proposed by Furlani. The magnetic gradient force is caused by an bias field and the 
polarized magnets and is found to be spatially uniform in the channel section which can be used for particle 
selecting or separation. We have derived the equations of nonlinear magnetization of magnetic particles 
which cause the harmonics of magnetophoresis. The Langevin model and generalized Clausius-Mossotti 
equation used show how the normal and longitude anomalous anisotropic effect the permeability of fer- 
rofluids, thus the magnetic force. Our analysis demonstrates the viability of using the microchannel system 
for various bioapplications and other characterization of fluid transporting and the time-varying magnetic 
field can be potentially used for an integrated magnetometer and influences the the viscosity and effective 
permeability in ferrofluids. 
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I. INTRODUCTION 



Nowadays magnetophoretic microsystems have been paid great attention in bio-technology for 
the integration of "micro-total-analysis"(/iTAS) UJ] because of its high degree of detection and 
selectivity. The biomaterials possessing low magnetic susceptibility can cause substantial contrast 
between the labeled and unlabled materials tM- Because of its polarization difference, the the 
particles exhibit rich fluid-dynamic behaviors such as magnetophoresis as well as various mag- 
netic responses. Magnetic cell separation can be applied using magnetic beads which coated with 
specific cell(core- shell microsphere), or the native magnetic susceptibility [13|,lZ|]. In some special 
cases such as blood cells, the red and white blood cells can be conducted using magnetophoretic 
separation based on their native magnetic properties: diamagnetic or paramagnetic Magnetic 
changes in red blood cells can also be used for separation of diseased cells ^5^. In such continuous 
microseparator, the ferromanetic wire(circular or square) put in close proximity under an external 
bias field cause strong magnetic field with high magnetic gradience. Miniaturized Cell separator 
can be integrated for various types of cell counting and collecting. The magnetophorsis with inte- 
grated soft magnetic elements have some advantages over electrophoresis with electromagnets [|6D 
for they consume no heat and cause no damage and other negative effect on the bio-cells. Furlani 
have recently demonstrated transport and capture behavior of magnetic particles such as Fe304 in 
the microsystem which consists of an array of integrated soft magnetic elements embedded under 
the microfluidic channel with also an external bias field. The elements can be polarized by the bias 
field, thus producing nonuniform field distribution which causes magnetophoretic force on mag- 
netic particles within microchannel. The cubic soft magnetic elements replace the wire, producing 
different separation and trap in geometry. 

In the present paper we emphasize on the characteristics of the magnetic composites(such as 
large magneticc constants and permanent moments) on the behavior of magnetophoresis in fer- 
rofluids and donot consider the equations for particle motion. In experiment the slow-moving 
transport is influenced by the viscous drag and thermal kinetics, thus the magnetic force can be 
measured in the quasi-equilibrium movement. 

Ferrofluids(or magnetic fluids) are colloidal suspensions containing single domain nanosized 
ferromagnetic particles dispersed in a carrier fluid [SJ. Since these particles can interact easily and 
form crystal-like structure in the presence of applied bias magnetic fields, which in turn can affect 
the viscosity and structural properties tremendously particles in ferrofluids have a wide variety 
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of potential biomedical applications such as label and manipulate biomaterials. The dynamic(ac) 
magnetic properties and magnetization-induced second-order harmonic generation are taken into 
consideration in the system. In experiment, the measurement of an ac complex magnetic sus- 
ceptibility of magnetic fluids is a suitable method to study the relaxation process of the magnetic 
dipoles of colloidal particles in magnetic fluids lioi 11]. The second-order harmonic generation is 
the phenomenon that the magnetization along the specific direction activates the originally silent 
tensor components for the second-order nonlinear optical susceptibility and is observed for sur- 
faces or interfaces of ferromagnetic materials(thin films) lll2ll . or polar antiferromagnets [|l3n and 
polar ferromagnets UM- 

The saturation of magnetic particles considered here will be different from Furlani's theory 
and Han et aVs experiment, which cause the nonlinearity to appear by two effects: normal sat- 
uration and anomalous saturation. In detail, the normal saturation arises from the higher orders 
of Langevin function at large filed region, and the anomalous saturation results from the particle 
chains with higher and lower dipole moments lisj] which is shifted under the influence of the 
field. The magnetic field inside the ferrofluid plays important role in the coupling between the two 
effects, which is similar in electrorheological and magnetorheological fluids. When suspension 
having nonlinear characteristic is subjected to ac magnetic field, the harmonics of magnetic sus- 
ceptibility can be induced to appear [16]. Our analysis demonstrates that the magnetic force will 
cause suitable separation across microchannel and will be affected by the anisotropic changes of 
magnetic dipole arrays. 
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II. THEORY AND FORMALISM 



A. Nonlinear magnetic moments in ferrofluids 

In the standard case the magnetic induction B is proportional with the field strength Hq, and 
have the relation B = /uHq, where /i is the linear permeability. At strong field intensities, non- 
linearities are introduced as B = /igHo + xHqUq where third-order and higher-order nonlin- 
ear coefficients are dropped and x and /ig are the nonlinear magnetic susceptibility and effective 
permeability for the longitudinal field case. Here we assume the nonlinearity is not strong and 
consider only the lowest-order nonlinearity for simplicity. 

In the ferrofluids, the average component in the direction of the field of the magnetic dipole 
moments can be expressed as (M^ ■ e), where (■ ■ ■ ) denotes (••■) = J j ^jyuie-^x/feflT dX\ e 
is the unit vector in the direction of the external field, and X stands for the set of position and 
orientation variables of all particles. Here u is the energy related to the dipoles in the sphere, 
and it consists of three parts: the energy of the dipoles in the external field Udc, the magnetostatic 
interaction energy of the dipoles Wmi, the non-magnetostatic interaction energy between the dipoles 
Unmi which is responsible for the short-range correlation between orientations and positions of the 
dipoles such as London- Van der Waals interaction energy. 

In this case, the effective permeability of ferrofluid /ig is determined by the generalized 
Clausius-Mossotti equation taking into consideration of dipolar interactions nisl lnll: 



At2 + fi'L(/ie - /U2) 3 V 3/cbT1 + i27r/r^ 
where ^2 represents the permeability of the host fluid, N the number density of the particles, ks 
the Boltzmann constant, T the absolute temperature, / the frequency of the applied magnetic field, 
r the relaxation time of the particles, and a the magnetizability of the particles. Here r can be ex- 
pressed as TbTn/ (Tf, + r„), where and r„ are the brownian relaxation time and the neel relaxation 
time respectively [|20|] . Our model describes the aggregation behavior in an external field by intro- 
ducing the longitudinal demagnetizing factor in Clausius-Mossotti equation, which deserves a 
thorough consideration: Eq. (2) should be expected to contain both u^i and Wnmi as (7^ is not equal 
to 1/3. For an isotropic array of magnetic dipoles, the demagnetizing factor will be diagonal with 
the diagonal element = 1/3. However, in an anisotropic array like ferrofluid, the demagnetiz- 
ing factor can still be diagonal, but it deviates from 1/3. In fact, the degree of anisotropy of the 
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system is just measured by how ql is deviated from 1/3. It is worth noting that ql < 1/3 in the 



present longitudinal field case. Furthermore, there is a sum rule for the factors, gL + 2gT = 1 nlsA . 



where gx denotes the transverse demagnetizing factor. Such factors were measured by means of 



computer simulations [|l9l]. Thus, to investigate the anisotropic structural information of the array, 
we have to modify the Clausius-Mossotti equation accordingly by including the demagnetizing 

factor. When we studied the field-induced structure transformation in ferrofluids, we can use the 

I > 

generalized Clausius-Mossotti equation [17] by introducing a local-field factor P' which reflects 
the particle-particle interaction between the particles in a lattice. 

The magnetic dipole moment nij satisfies Langevin function = ms(coth7 — ^) where 
7 = moHo/ (ksT) and denotes the saturation magnetization of particles. For the whole array 
of magnetic moments, = Xir=i(™*(/^oo + 2/i2)/3/i2)i where fioo represents the permeability 
at frequencies at which the permanent dipoles cannot follow the changes of the field but where the 



atomic and the electronic magnetization are still the same as in the static field [|15l1 . Therefore, fi^o 
is the permeability characteristic for the induced magnetization. In practice, jjoo can be expressed 
in the expression containing an intrinsic dispersion, 

/.oo = /i^(O) + (2) 

where /ioo(O) is the high-frequency limit permeability, and Afi stands for the magnetic disper- 
sion strength with a characteristic frequency fc- Harmonics of magnetic moments can be obtained 
through Frohlich model bill by 



M..e^-J^. (3) 

Here Hp gives the magnetic field inside the spherical situated in medium with permeability fXe, 
and has the form — Hq + 3x^°° so HL Thus taking into account the higher derivatives of the 
average moment we obtain 



Noticing the expression for (M^^ ■ e) and Eq. (3), it is easily derived 



^ -(M,.e)U,=o = r^[((M,.e)2)-(M,-e)2]U,=o = r^(Mj)o, (5) 
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+6((M,-e)2)(M,-e)2-3((M,.e)y 
-6(M,-e)^ + 6(M,-e)2((Md-e)2) 

-{M,-e){{M,-e)')]\H,=o 

[3(Mi)o-5(M,2)^]. 



[((M,-er)-3((M,-e)3)(M,-e) 



15(A;bT)3 
Using the same method, we obtain 



(6) 



(Mi)o 



3/i2 



ni 



y V 3/i2 



ni 



ni 711 



[yPo ^{cos Oij J2 cos Ors)]- 



(7) 
(8) 



r=l s=l 



For the numerical calculation, we have 

^(cos%) = 1, 

"2 "2 "2 2 

^(cos% ^ ^cos6'r5) = -(5^1-2) 
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(9) 
(10) 



j=l r=l s=l 

In the presence of external oscillating time-varying magnetic field cite [|22l1 . the magnetic parti- 
cles will have nonlinear characteristics. In the experiment, the second term or higher order terms 
of magnetization or force can be obtained using mixed-frequency measurements [123'] . When we 
apply an external field such as i/o(^) = -f^dc + Hs.c{t) = H^c + i/acsincut, the orientational 
magnetization Mo will contain harmonics as 



Mo = M^<^') + sin ut + M2<^ cos 2ujt + M^^ sin 3ujt + ■ ■ ■ . 



(11) 



Here cu = 2nf and Hdc denotes the dc field which induces the anisotropic structure in the 
ferrofluids, and H.^c{t) stands for a sinusoidal ac field. Applying Eq. (11) into Eq. (4), (5) and (6), 
after tedious calculation, the harmonics of magnetic moments can be expressed as 
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^) = ifdc Ji + IhIH,,x + Hlx, (12) 

= H^.h + \hIx + ^H^cHlx, (13) 

= -^i/^^i/dcX, (14) 

Ms. = -\hIx, (15) 

with Ji = /ie — /Xoo- Below we set ifac = -f^dc = -f^ for simplicity, thus all harmonics for Mq can 
all be expressed as Mq = x H- 



B. Magnetophoresis in microchannel 

Now we will investigate the magnetophoretic behavior in microchannel. In a standard case, 
a magnetically polarizable object will be trapped in a region of a focused magnetic field, pro- 
vided there is sufficient magnetic response to overcome thermal energy and the magnetophoretic 
force [24I. For a magnetically linear particle under magnetophoresis, the effective magnetic dipole 
moment vector induced inside takes a form very similar to that for the effective dipole moment 
of dielectric paricle, m = 47ri?|^gMoRe[CMF], where i?cfr is the effective radius of spherical or 
spheroidal particle and CMP is the Clausius-Mossotti factor along the direction of external field. 
It is noted that particles are attracted to magnetic field intensity maxima when 6 > as positive 
magnetophoresis and negative magnetophoresis corresponds to 6 < 0. For spheroidal particle with 
permeability fi2, CMF = | ^ +g^^(~^^2At ) • magnetophoretic force exerted on the particle in a 
nonuniform magnetic field H can be written as, 

i4 = 47r/iei?^ffRe[CMF]Mo ■ VH (16) 

The microsystem consists of one integrated soft-magnetic elements which is embedded in a 
nonmagnetic substrate beneath a microfluidic channel as shown in Fig. (1), in which the magnet 
is 2w wide and 2h high. The magnetic particles in ferrofluids which pass through the channel can 
separated according to the different field gradient distribution, and the nonmagnetic particles will 
be rinsed away. Y axe corresponds to the longitude case mentioned above. Below we will show 
that when the number of magnetic element increases, the magnetophoresis will become more 
uniform in both parallel and perpendicular direction to that of particle transport. The bias field 
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-^bias along y axe and field from the magnet i?mag which is saturated(saturated magnetization Mes) 
when the bias field is added, both contribute to the magnetic field exerted on the particle. Under 
high frequency of bias field(> lO'^Hz), the magnetization of magnet can be viewed constant. Thus 
the magnetophoretic force is decomposed into two components: 

F^^ = 47r/i,i?,\Re[CMF]x'[i^n.ag,x^^^ + {H^.,,y + i7bias,y)^^^], (17) 
F^, = 47r/iei?e'frRe[CMF]x'[iJ„.ag,x^^^ + (^mag,y + i^bias.y) ^^^] • (18) 

where H^^g can be expressed as iY^ag.x = ^{H ^^XZy^t^S ^ " ^"^^WW+StW^^ 
i^mag,y = ^ {tan"^ [ (^^^^l^+yJ.f,^ ] " tan'^ [ ^J^l2~y2lh^ ] } from fel]. The magnetophoretic force 
across the channel can then be be accurately calculated. It could be noted from the equations above 
that without the bias field magnet-element cannot be magnetized and bias field contributes to the 
magnetophoretic force, but the bias field itself produce no field gradient in the microchannel. 
Because the importance of the the contribution from the bias field which causes harmonics of 
magnetization in the particle, below we will calculate it alone. 



III. RESULTS 



Now we are in the position to study the behavior of magnetophoresis in ferrofluids. Specif- 
ically, the external field is set to be H^c = H^c = O.OlA/m, high-frequency limit permeabil- 
ity //oo(0) = Stt X lO^^H/m, dispersion strength A/i = 32n x 10^''H/m, medium permeability 
fX2 = 4:71 X lO^^H/m, permanent magnetic moment for each particle tuq = 10^^'^A/m, characteris- 
tic frequency fc = 5x lO^Hz, the relaxation time r = 5 x 10^ ''s, the number density = 10^m~^, 
the length of magnet w = h = 50/im, Boltzmann constant ks = 1.38 x 10~^'^J/K, the effective 
radius i?cfT = 500nm and the average magnetization of single particle a = lO^^m^'^. 

Figure 2 displays the temperature and anisotropic effect on CMF in magnetophoresis. In Fig. 
2(a), it is shown that increasing the temperature T causes slight increase in CMF(about 3%), thus 
has little effect in magnetophoresis. From Fig. 2(b), the magneophoretic force exhibits strong 
sensitivity to the anisotropic factor g^. When CMF is negative(particles will be repelled from 
the maximum magnetic gradient), the larger is, the smaller repulsion will be. When CMF is 
positive, with increaseing gi such an attractive force becomes stronger. Furthmore, we predict 
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the anisotropic-dependent crossover frequencies at which there is no net force on the cell particle. 
The crossover frequency is monotonically increasing function of ql, dependent on whether the 
variation of magnetophoretic force is negative or positive. 

Figure 3 and 4 shows the oriental, fundamental, second- and third- order harmonics of the 
magnetic force and Fy induced by the bias field mentioned above as a function of the field 
frequency for various g'L in the longitudinal field case. The particle is put in the position x=50nm 
and y=300nm(see the coordinate in the inset of Fig. 5). The harmonics of magetization of particles 
will change accordingly as the system alters from isotropic case {ql — 1/3) to anisotropic {ql ^ 
1/3) because of the appearance of the particle chains. In detail, stronger anisotropy (namely, 
decreasing the longitudinal demagnizing factor gi) leads to larger magnetic force in the low- 
frequency region. It can also be observed there is two plateaus in the range 0.03-3MHz and 10^- 
lO^MHz where the magnetic force have slight change and the second and third order harmonics 
of magnetic force are negligible compared with lower order ones. In fact for the transverse field 
case, it could be concluded from the sum rule between ql and qt that ql + '^^Qt — 1- Because 
of the coupling between the applied dc and ac magnetic fields, the even-order harmonics are also 
induced to appear besides the odd-order harmonics for the longitudinal field case, even though 
only the cubic nonlinearity is considered due to the virtue of symmetry of the system. In addition, 
the harmonics shown in Figs. 3 and 4 are nonzero dX ql — 1/3 is because of the presence of 
external fields (i.e., u^^ 0), even though there is no particle interactions (i.e., u^i — linmi = 0) 
as g'L = 1/3, the nonlinear behavior due to the normal saturation could still be induced to occur. 

Fig. 5 and 6 displays the magnetic force which is spatially varying across the channel section, 
when one and two magnets are embedded in the substrate under frequency lOKHz of magnetic 
field. In Fig. 6 the two magnets are embedded 250/im away and the x — y coordinate is placed 
in the midpoint of magnets. Comparing Fig. 5(A) with Fig. 6(A), we find that the force F^^ 
for the purpose of separating magnetic particles which is attractive in some regions and repulsive 
in others along the x axe becomes more uniform, and the force Fy has similar nature. In detail 
the absolute value of magnetic force F^ is symmetric in the x direction and become stronger on 
the edge of the microchannel, while Fy is always negative which will repel the particles from 
the magnets. The uniform separation technique will cause the magnetic particles with different 
size or permeability become apart into layers in geometry which can in fact not only be used in 
ferrofluids. The calculation also shows that when more magnets are uniformly embedded, the 
separation of particles are more efficient. It should be pointed out that the distance between the 



9 



nearby particles is important to cause spatially uniform magnetic force and it can be chosen for 
different purpose of separation which may be determined by the capture efficiency for a specific 
particle sorting or transporting. The time-varying magnetic field can be potentially used for an 
integrated magnetometer and influences the the viscosity and effective permeability in ferrofluids. 

IV. CONCLUSION 

We have studied the magnetophoretic particle separation of ferrofluids in microchannel and 
its nonlinear behavior. The magnetic gradient force is caused by an bias field and the polarized 
magnets and is found to be spatially uniform in the channel section which can be used for particle 
selecting or separation. We have derived the equations of nonlinear magnetization of magnetic 
particles which cause the harmonics of magnetophoresis. The Langevin model and generalized 
Clausius-Mossotti equation used show how the normal and longitude anomalous anisotropic effect 
the permeability of ferrofluids, thus the magnetic force. Our analysis demonstrates the viability 
of using the microchannel system for various bioapplications and other characterization of fluid 
transporting. 
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Figure Captions 

Fig. 1. (Color online) Schematic graph showing one integrated soft-magnetic elements embeds 
in a nonmagnetic substrate beneath a microfluidic channel through which ferrofluids flow. 

Fig. 2. (a) CMF for different temperatures vs the frequency f of magnetic fields, (b) CMF vs 
the frequency fof magnetic fields for different anisotropic factors gi- 

Fig. 3. Oriental, Fundamental, second and third order harmonics of the magnetic force vs 
the field frequency for various g'L in the longitudinal field case. 

Fig. 4. Same as Fig. 3, but for another force component Fy. 

Fig. 5. (Color online) (A) The spatial magnetic force F^ across the channel section, when 
one magnet embeds in the substrate under frequency lOKHz of magnetic field. (B) The spatial 
magnetic force Fy across the channel section. 

Fig. 6. (Color online) Same as Fig. 5, but for two magnets case. 
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FIG. 1 : Jian, Zhang and Huang 
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FIG. 2: Jian, Zhang and Huang 
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FIG. 3: Jian, Zhang and Huang 
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FIG. 4: Jian, Zhang and Huang 
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FIG. 5: Jian> Zhang aiid Huang 
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